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ABSTRACT 

The low frequency radio emission of starburst galaxies is informative, but it can be 
absorbed in several ways. Most importantly, starburst galaxies are home to many H 
II regions, whose free-free absorption can obscure low frequency radio waves. These H 
II regions are discrete objects, but most multiwavelength models of starbursts assume 
a uniform medium of ionized gas, if they include the absorption at all. I calculate the 
effective absorption coefficient of H II regions in starbursts, which is ultimately a cross 
section times the density of H II regions. The cross section can be easily calculated 
by assuming that H II regions are Stromgren spheres. The coefficient asymptotes to a 
constant value at low frequencies, because H II regions partially cover the starburst, 
and are buried part way into the starburst 's synchrotron emitting material. Consid- 
ering Stromgren spheres around both O stars and Super Star Clusters, I apply the 
calculations to the low frequency radio spectrum of M82. Far from being opaque at 
low frequencies, I find that M82 mostly transmits its radio flux. I also find that star- 
bursts are transparent down to a few MHz to other possible absorption processes, such 
as free-free absorption from the diffuse superwind phase, synchrotron self-absorption, 
and the Razin effect. Hence, starburst galaxies should be observable with new low 
frequency radio telescopes. 

Key words: radio continuum: general - radio continuum: ISM - galaxies: starburst 
- H II regions - galaxies: individual (M82) 



1 INTRODUCTION 

The radio spectrum of star-forming galaxies, including star- 
bursts, is dominated by synchrotron emission from cosmic 
ray ( CR) electrons and positrons (e ) in diffuse magnetic 
fields |Condo3[l99i). There is also free-free emission, which 
typically comes from H II regions in a galaxy. Both of these 
emission processes are associated with star-formation: cos- 
mic rays are generated somehow by star formation (pos- 
sibly through shock acceleration in supernova remnants), 
and H II regions surround young, massive stars that pro- 
duce ionizing radiation. Since synchrotron emission has a 
steeply falling spectrum (typically Sl ynch oc v~ ' 7 ) whereas 
free-free emission does not (Sv oc v~ 0A ), th e synchrotron 
emission dominates below about 30 GHz (e.g.. lCondon|[l993 : 
iNiklas. Klein. Wielebinskilll997l h 

There is a wealth of information available at MHz 
frequencies. In particular, different cooling processes may 
set the CR e lifetime at different energies. At low fre- 
quencies, bremsstrahlung, with an energy-independent 
loss time, and ionization, which is most effective at low 
energies, become more important than synchrotron and 



Inverse Compton coolin g , which grow stronger at high 
energ ies (iHummell Il99ll; [Thompson et all 120061 ; iMurphvl 
12009 : lLacki. Thompson, fc Quataertl [2010). Any escape, 
whether diffusive or advective, also becomes more impor- 
tant relative to synchrotron at low frequencies. Therefore, 
the synchrotron radio spectrum should flatten at low 
frequencies and steepen at high frequencies. Indeed, this 
behaviour is often seen in the radio spectra of starburst 



galaxies (e.g ., Clemens et al.ll200 i lWiUiams fc BowerllioTol : 
iLerov et al.l 12011 ). The detailed radio spectra is then 
useful in constructing models of the cosmic ray population, 
helping to con strain the poorly understood magnetic field 



stren gth (e.g.. iTorres! 2004 : Domingo-Santamarfa fe Torres! 
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sion might also reveal the presence of a 'pion bump' in the 
e spectrum, in which the spectrum of secondary e from 
pion decay falls off due to t he kinematics of pi on production 
in proton-proton collisions (|Rengaraianl 2005). 

In practice, our ability to understand the low frequency 
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radio spectrum of starbursts is limited by other processes 
which alter the radio spectrum. The most important is 
free-free absorption by ionized gas in galaxies. The Galac- 
tic radio spectrum has a turnover at ~ 3 MHz, largely 
caused b y free-free absorption i n the diffuse Warm Ionized 
Medium (lAlexander et~atlll969l ; iFleishman fc Tokarevll 19951 ; 
iPeterson fc Webberl20ol ). Free-free absorption is even more 
important along lines of sight through dense H II regions, 
whic h can become optically thick even at GHz frequencies 
(e.g., iMcDonald et al.ll2002f) . Free-free absorption, like ion- 
ization and bremsstrahlung cooling, can flatten the low fre- 
quency radio spectra of starbursts. The spectral curvature of 
starburst galaxies has therefore also been interpreted as free- 
free absorption (|Klein. Wielebinski. fc Morsil Il988l ; [C 



iii 



ll99d : IClemens et al.ll20ir3 T 

The observational prospects for low frequency studies 
of star-forming galaxies are good. There is increasing inter- 
est in low frequency radio observations, due to their value 
in observing high redshift 21 cm lines, among other rea- 
sons. The Giant Metrewave Radio Telescope (GMRT) is 
specifically designed to provide interferometric data for ra- 
dio sources in the 50 MHz to 1.5 GHz range with high 
sensitivitjQ. The 74 and 333 MHz systems on the Very 
Large Array (VLA) provide high angular resolution images 
(Kass ini et al.l l2007t). and have completed a survey of the 
northern sky (|Cohen et al.l [2007b . It is currently being up- 
graded for use on the Jansky VLA. The Low Frequency Ar- 
ray (LOFAR), now coming online, is a new radio telescope 
with long baseline interferometry capabilities, and can go all 
the way down to 15 MHiQ The Leiden LOFAR Sky Surveys 
Project will image some nearby star-forming galaxies in ra- 
dio, possibly including starbursts like M82. LOFAR will be 
joined by the 21 cm pathfinder experiments at frequencies 
above 100 MHz and possibly the Long Wavelength Arr ay at 
frequencies of 10 to 88 MHz fe.g., lEllingson et al.ll2009h . Ul- 
timately, the Square Kilometre Array (SKA) should be able 
to observe galaxies down to 70 MHz with high sensitivitjQ. 

However, there has been relatively little work on the 
theory of the low frequency radio spectra of starburst galax- 
ies. For model fitting, the typical assumption, if free-free 
absorption is even considered at all, is the uniform slab 
model, in which both the free-free absorption and emis- 
sion comes from a uniform density ionized medium pervad- 
ing the synchrotron-emitting region (for examples of uni- 
form slabs used to fit s tarburst radio spec t ra, see, e.g. , 
Sopp fc Alexander! Il99ll; ICondon et all Il99ll; ICariih] Il99rj; 



Torre s 2004; Clemens et al.ll2010l ; William s fc Bowerl 



2010h . 



The best measurements are at GHz frequencies, where the 
integrated free-free absorption is often expected to be small 
and the details of the absorption may not matter much, 
though it has been claimed to be important at GHz frequen- 
cies in Arp 2 20 and other Ultraluminous Infrared Galax- 
ies (ULIRGs) (ICondon et al.lll99ll ; ISopp fc Alexanderlll99ll ; 



IClemens et alj|2008l l2010h . However, from a theoretical 
point of view, this approximation is likely to be too sim- 
ple: in the Milky Way, while the free-free absorption largely 



1 See http://www.gmrt.ncra.tifr.res.in 

2 At http://www.lofar.org. 

3 With a home page at http://lofar.strw.leidenuniv.nl/. 

4 http://www.skatelescope.org. 



comes from the diffuse ionized medium, the free-free emis- 
sion comes from compact H II regions. On the other hand, 
in starburst galaxies, at least part of the free-free absorption 
is observed to co me from discrete H II regions, as seen in the 
Galactic Centre dBrpgan et al.|[2003l : lNord et al.ll2006t ), and 
in M82 l| Wills et al.l Il997f) . Discrete H II regions partially 
cover starbursts: some of the sightlines through the star- 
burst do not intercept an H II region and are unabsorbed at 
any frequency. 

The alternative is to consider starbursts as a collection 
of discrete H II regions. A simple version of this approach 
has been considered in the context of radio recombination 
line studies, where the radio spectrum is needed to calcu- 
late stimulated emission (e.g.. [Anantharamaiah et al.l 1 19931 : 
IZhao et all Il99rj ; iRodriguez-Rico et alj|2005l ). The typical 
assumption in these models is that all of the H II regions 
have the same density, temperature, and radius, and are 
located in the midplane of a starburst disc. Each H II 
region then shadows the synchrotron-emitting region be- 
hind it: with this assumption, the free-free abso rption can 
be predicted (e.g.. lAnantharamaiah et al.| [l993). However, 
the model does not work well when the H II regions are 
both optically thick and have a covering fraction near 100% 
|Anantharamaiah et al.| [l993). and it does not allow for H 
II regions of different radii. 

In this paper, I study effects on the low frequency spec- 
tra of starbursts. I focus on free- free absorption, the most im- 
portant process. I calculate the amount of free-free absorp- 
tion by assuming it comes from discrete H II regions around 
stellar clusters in starbursts. In essence, I treat a starburst 
as a pudding with Stromg ren spheres mixed in. My mo del 
generalizes the approach of lAnantharamaiah et all (|l993l ) to 
allow H II regions of different radii, here assumed to corre- 
spond to young star clusters of different masses, and allows 
for optically thick H II regions with high covering fraction. 
I make spectral fits using these model to the low frequency 
radio spectrum of M82 in Section [4] accounting for the pos- 
sibility of spectral curvature. I also consider whether other 
processes, such as free-free absorption from the keV plasma 
of the starburst wind, the Razin effect, or synchrotron self- 
absorption, could cut off the low frequency spectrum in Sec- 
tion O 



2 A METHOD FOR CALCULATING 
FREE-FREE ABSORPTION FROM 
DISCRETE H II REGIONS 

The nonthermal continuum radio emission of starbursts are 
thought to pervade the entire starburst region, because the 
radia ting cosmic rays can diffuse from their acceleration sites 
(e.g., iTorres et al1[2012h . The H II regions are intermixed 
in the starburst region, so in this sense the uniform slab 
model is correct. However, a model of truly uniform ion- 
ized gas assumes the free-free absorption comes from a low 
density, high filling factor medium, which would have a low 
turnover frequency, but would be highly opaque below the 
turnover frequency. However, a more appropriate assump- 
tion is that there is a uniform density of H II regions rather 
than of ionized gas. A uniform collection of discrete H II 
regions is a high density, low filling factor medium, with a 
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Figure 1. H II regions (rilled circles) in a slice of the starburst 
with area A and thickness As absorb radiation from the column 
behind it. The column can be small relative to the rest of the 
starburst (shaded in grey), but each slice is large enough to con- 
tain a representative sample of H II regions. A beam of light is 
emitted from the background column, enters perpendicular to the 
slice with intensity Z? n , and emerges with intensity Z£"' k , where 
I is flux per solid angle. The fraction of light absorbed passing 
through the slice, which can be calculated by summing over in- 
dividual H II regions in the slice, is directly related to the mean 
effective absorption coefficient of the starburst from H II regions. 
The coefficient then can be used in a uniform slab model (as sug- 
gested by the uniform grey shading of the rest of the starburst). 



higher turnover frequency, but also translucent below that 
frequency. 



2.1 Derivation of Absorption Coefficient for 
Discrete H II Regions 

I start by considering how a thin slice of the starburst with 
thickness As absorbs background radio emission passing 
perpendicularly through it (Figure [l}. The slice has an area 
A and thickness large enough to contain a representative 
sample of the starburst's H II regions, but otherwise can be 
small compared to the starburst. I ignore the H II regions' 
effects on radio emission emitted within the slice - as long 
as the slice is thin, and as long as the H II regions have a 
small filling factor, this should be a valid approximation. If 
the filling factor is large, a true uniform slab model is more 
accurate anyway. 

Suppose that there are M H II regions (indexed by the 
number m). Each H II region passes a fraction <I> m of the 
normal incident radiation I^l c ^A on the slice, where I is flux 
per solid angle. For example, if a totally opaque sphere of 
radius r sat in a slice of area A, $ = 1 - (nr 2 /A). Then the 
background radio intensity after passing through the slice is 



n 



<A. 



(1) 



Taking the logarithm of both sides, we can rephrase this as 
an optical depth: 



M 

£ 

m— 1 



In 



In 



(2) 



Because each H II fraction only covers a small portion of 
the starburst, a convenient approximation is to define 4> m = 
1 — and then to assume that \n§ m ~ — <f) m . This is 
equivalent to saying that the H II regions do not overlap 
within the slice. Then we have 



M 
m — 1 



Ar n , 



(3) 



If there are J types of H II regions, where each type of 
H II region absorbs the same amount of background radio, 
and if there are Kj H II regions of each type j within the 
slice, then 



Ar ci 



J 

£ 

3=1 



Ka 



(4) 



More realistically, instead of having a few distinct types of H 
II regions, there will be a distribution function parametriz- 
ing the absorption properties of H II regions. Suppose the 
number of H II regions with some relevant physical quantity 
(such as radius) between q and q + dq and within a vol- 
ume dV is dN j (dqdV) . Then, since the volume of the slice 
is AAs, the effective optical depth is 



Ar c; 



dN 
dqdV 



(AAs)<j>(q)dq. 



(•») 



Of course, this can easily be generalized to more than one 
parameter. 

The effective optical depth across the slice translates to 
an effective absorption coefficient a c s = Ar e ff/(As) of the 
starburst. Furthermore, the quantity A<f>(q) is an effective 
absorption cross section a for each H II region. So we have 



Qoff 



dN 
dqdV 



<r{q)dq. 



(6) 



Therefore, the absorption coefficient reduces to a number 
density times a cross section, as might be expected. 

If the density of H II regions is constant throughout 
the starburst, the equation of radiative transfer through the 
starburst, 



dl_ 

ds 



-a e «I + j 



has the uniform slab solution: 



3 

Otett 



(7) 



(8) 



except that the optical depth is the effective optical depth 
defined as r e a = OeffS, where s is the sightline length through 
the starburst. Of course, the density of H II regions may 
itself vary, for example, decreasing towards the edge of the 
starburst. However, the basic principle remains the same: 
the effective absorption coefficient is calculated at each point 
using equation [6l and then used in the radiative transfer 
equation to find the synchrotron intensity on each sightline. 
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2.1.1 The effective cross section for translucent spheres 

I treat H II regions as Stromgren spheres, which free-free 
absorb the radiation behind them. The effective cross section 
of a sphere of radius Rs and absorption coefficient ami is 
equal to the projected area of the H II region times the 
fraction of background flux it obscures: 

rRs r 

l_ e - 2a H"V fl s 



2-rvy 



■nR%- 



2n- 2 

ZQ HII 



ti- 



dy (9) 
- 2 " H " Rs (l + 2oHiiiJ S )] (fO) 



When ami is very large (at low frequency, for example), 
the spheres become totally opaque and the effective cross 
section is just irR s . If instead ami is very small, as at high 
frequency, the effective cross section goes as (4/3)7rami^?s- 
A simpler approximation, the blue edge approximation, 
for the effective absorption cross section of the sphere as- 
sumes it is completely transparent when thii ~ amiRs < 1 
and completely opaque when thii > 1: 



irRs {amiRs > 1) 







(amiRs < 1) 



(11) 



This approximation works well at high and low frequencies, 
but at intermediate frequencies, it can fail. 

2.1.2 The covering fraction 

The calculation of the covering fraction of the H II regions 
is analogous to the calculation of the absorption of the radio 
flux. Essentially, we are interested in how much of the back- 
ground sky is blocked by the H II regions. We can phrase 
this in terms of a flux calculation: if H II regions are totally 
opaque, the covering fraction is equal to the fraction of a 
uniform background intensity they absorb. Thus, we have 
for a distribution function of H II regions 

For the covering fraction calculation, the H II regions are 
treated as a foreground screen rather than as a uniform slab, 
as it does not matter how deep into the starburst any H II 
region is. From the solution to absorption from a foreground 
screen, the probability any given sightline is covered is 

CKcovcrS (13) 



1 



where s is the length of the sightline, assuming a covor is 
uniform throughout the starburst. 

The length of a sightline may vary; for example, in an 
edge-on disc, it is shorter towards the rim than through the 
centre. The covering fraction for the entire starburst is the 
average of the covering probability P CO vcr over all sightlines 
that pass through the starburst: 



=J-f 

^obs J 



(i 



')dfio 



(14) 



In this equation, Slobs is the solid angle subtended by the 
entire starburst from Earth. Note that since starbursts are 
small on the sky, Q b B is directly proportional to the pro- 
jected area of the starburst A pTO j : 



/cover — ~Z / (1 e covor )dAp 



(15) 



A natural consequence of / CO vcr < 1 is that we expect 
the optical depth to increase to some maximum value at 
low frequency and no further. If there are just a few H II 
regions, then the covering fraction is small and parts of the 
starburst remain completely unobscured. However, even if 
the density of H II regions is high, discrete H II regions 
tend to overlap from our point of view and cover each other 
instead of background material. This leads to a covering 
fraction less than 1. Thus, unlike a diffuse ionized medium, 
discrete H II regions lead to non-zero radio emission at low 
frequencies. 



2.1.3 The filling fraction 

Likewise, we can calculate the filling fraction from a distribu- 
tion of Stromgren spheres. The starburst has a volume Vsb, 
and a Stromgren sphere of radius Rs has volume 4/37rJ?|. 
Therefore, each sphere leaves a fraction 1 — AirR s /(3V) of 
the starburst unfilled, and the filling factor can be calculated 
using 



1 - /flu « exp 



dN ( 4nR s (q) 3 



dq \ 3Vsb 



dq 



(16) 



When the filling factor is small, the spheres do not overlap, 
and 



/flu 



f dN_ ( ^Rs{qf 
J dq \ Wsb 



(17) 



Since Rs is much smaller than the size of the starburst, 
the filling factor of H II regions can be much smaller than 
1 even if the covering fraction is nearly 1. As a result, a 
sightline through the starburst will not immediately inter- 
cept an H II region, even if the covering fraction is nearly 
1. Instead, the first absorbing H II region on the sightline 
is buried past unobscured synchrotron-emitting space: this 
ensures that there is still synchrotron radio flux at low fre- 
quencies even if the covering fraction is high. 



2.2 Using the absorption to calculate integrated 
flux 

The uniform slab solution provides the intensity along a sin- 
gle sightline, which can be appropriate if the starburst is 
resolved. However, we often are interested in the integrated 
flux of the starburst. The sightline will have different lengths 
as it passes through different parts of the starburst. I derive 
for the reader how much integrated flux remains after free- 
free absorption for three common geometries: a disc viewed 
face-on, a disc viewed edge-on, and a sphere. 



2.2.1 A face-on disc 

The flux observed at Earth is proportional to the luminosity 
emitted by the starburst into a unit solid angle, dL/d£l crn = 
J IdAproj . For a face-on disc of radius Rsb and midplane-to- 
edge height hsB (edge-to-edge height 2/isb), this is simply 
dL/dQ C m = TrRg B I, or 



dh 



(1- 



-2a e ff hSB 



(18) 
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If the starburst had no absorption, it would have dL/d£l cra — 
2TvRg B hsBj- Therefore, the ratio of actual flux to flux with- 
out absorption is 

if— = ~^r- h • ( 19 ) 

-funabs ^«Gff«SB 

asymptoting to 1 — a e fr/isB when a e ff is very small and the 
starburst is nearly transparent. 

The covering fraction for this geometry is just 

_ y _ e -2a CO v(!rhsB 



2.2.2 An edge-on disc 



(20) 



If the starburst disc is instead observed edge-on, the sightline 
length s through the disc varies between and -Rsb. Defining 
y as the impact parameter from the centre of the disc, so 
that y = yj i?g B — s 2 /4, the emergent flux is related to 



dL 



2h SB x 2 



3 

OeB 



(I - e-"»« s ) dy (21) 



Ahsi 



tteff 



Rse 



The ratio of actual flux to flux without absorption is 



F 



-Frrrrabs 7TQ c ffi?SB 



-2a e ffHsB\A : 



du 



(23) 



When the disc is nearly transparent, with small a e s, the 
ratio is approximately 1 — 8a or r.RsB/(37r). 

The covering fraction is 1/(4/isb-Rsb) f*! B f' l f B (1 - 

' x ' J — USB J —"SB 



e acova ' s )dxdy, which simplifies to 



-2a CO vorflsB\/ 1-« 2 



/cover = 1 - / e-* a ™*"* nSB V — du. (24) 

Jo 

When a COVC r-RsB < 1, /cover « (7r/2)i? S B« cover • 



2.2.3 A sphere 

Suppose instead a starburst is a sphere with radius 7?sb- 
Much like the case of an edge-on disc, the sightline length s 
varies between and J?sb , depending on the impact param- 
eter y from the centre of the sphere, with y = y/ R^ B — s 2 /4. 
The emergent flux is now proportional to 

= 2tt f SB y^-(l-e- a " aa )dy (25) 
aeff v 



dL 



dQ c 



CtefC 



Rsb + 



( Rsb 
\ aeff 



+ 



2a 2 fr 

err 



-2a off R SE 



2a 



§96 1 



In a transparent sphere, the luminosity per solid angle is 
dL/dQ cln = (4/3)7ri?g B j. Thus the ratio of actual flux to 
unabsorbed flux is 



F 



Fun a 



4a c ff-RsB 



,2 , (Rsb . 1 



#SB + I — + _ , 

V aeff 2a^ fl 



-2a c ff-RsB 



1 



2a 2 ff 

err 



In the limit when a e ff is small, this ratio is 1 — (3/4)a e ff-RsB- 
The equation for the covering fraction is now / CO ver = 
l/(ir.R§ B ) f* SB 2-Ky{l-e~ a ^ a )dy. This can be solved, giv- 
ing: 

1 r -2r»_„, 



/ 



1+ 



2R* 



'(2a C over-RsB+l)-l].(28) 



This reduces to f c , 



(4/3)a CO ver-R when a C over-RsB -C 1. 



3 APPLICATION TO POPULATIONS OF H II 
REGIONS IN STARBURSTS 

I now calculate the free-free absorbing properties of H II 
regions in starbursts. I assume that the ionizing sources 
in starbursts are either O stars or stellar clusters acting 
as point sources, each with an ionizing photon luminosity 
Qi on . I then assume the H II regions are spheres with the 
Stromgren radius: 



Rs = 



( 3Q io 



1/3 



(29) 



The recombination constant oib is equal t o a_g = 2.56 x 
10 _13 (T/10 4 K)-°' 83 cm 3 s- 1 (|Draindl20lfal ,. 

For a lon g -lived (^10 Myr) continuous starburst, 
iLeitherer etafl i| 19991 ) finds the ionizing photon luminosity 



Qio 



2.18 x 10 53 ph s" 1 



SFR 



M 



(30) 



assuming Solar metallicity and a Salpeter mass function 
from 0.1 to 100 Mq. I divide this ionizing luminosity into 
a calculated number of O stars or star clusters to find the 
typical ionizing luminosity of a source. 



3.1 H II regions around individual O stars 

One possible assumption is that the ionizing sources are in- 
dividual O stars. A typical O star has an ionizing photon 
luminosity of ~ 10 49 ph s _1 . Given the total ionizing pho- 
ton flux of eqn. 1301 the number of ionizing O stars in the 
galaxy is 



No = 21800 



SFR 



Qto 



Mq yr- 1 J yiO 49 ph s" 
The Stromgren radius around each O star is 



-Rs = O.i 



pc 



10 49 ph s" 1 



1/3 



nil 



1000 cm- 



-2/3 



(31) 



(32) 



for T e = 10 4 A". The free- free of a fully ionized hydrogen 
plasma (with n e = uh) is 



ami = 0.018 cm" 1 



niT 



- 3/2 ^Vff 



cm- 6 K- 3/2 Hz- 



(33) 



with the Gaunt factor git usually havi ng a value near ~ 10 
for 1 4 K plasma at MHz frequencies (|Rvbicki fc Lightmanl 
ll979T ). Therefore, each H II region becomes optically thick 
(ohii-Rs = 1) when 



, T v -3/4 . s 2/3 

m = 615 MHz (_) (■ m ^) 

1/2 



10 49 ph s- 



1/6 



(34) 



If the Stromgren spheres are in a uniform density 
medium, then all of the spheres are the same size, and it 
is straightforward to calculate the effective absorption coef- 
ficient for a starburst disc: 
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Figure 2. Effective absorption coefficient for H II regions of var- 
ious densities, as normalized to a volumetric star-formation rate 
of 1 Mq yr — 1 kpc -3 . This plot shows the case when each O star 
is single and has a Stromgren sphere around it. 



«efl = 



No 



2tt-RH b /isb 



(35) 



In the blue edge approximation, the result is simple: 

( Rsb 



Ctctt 



(1.2 kpc) -1 



Qto 



SFR 



M yr" 1 J \250 pc 



hsB 
100 pc 



10 49 s" 1 



-1/3 



ii ii 



1000 cm- 3 



-4/3 



(36) 



when v<vs- For star- formation rates of 10 Mq yr -1 , similar 
to those found in M82, we can already see that there will 
typically be several hundred parsecs of synchrotron emitting 
material on any line of sight before hitting an H II region. 
Since the typical line of sight through M82 has length ~ 
Rsb ~ 250 pc, there is about one Stromgren sphere on any 
line of sight, if most of the massive stars are isolated. 

In Figure [2] I show the effective absorption coefficient 
from Stromgren spheres around O stars for densities 1 — 
10 6 cm -3 . The absorption coefficients reach a plateau at low 
frequency, with a value that depends on hydrogen density. 
At high frequencies, the absorption coefficients for different 
hydrogen densities all have the same value. In this case, the 
cross section of Stromgren spheres with volume Vs is Vsami 
(section Pi .1.1 p . which does not depend on density. 



3.2 H II 



regions around Super Star Clusters: A 
simple model 



Much of the star-formation in starb urst galaxies occurs in 
bound super star clusters (SSCs; e. g., lO'Connell et al.lFl995l : 
iMelo et al.ll2005l ; ISmith et alj|2006ft . The SSCs have a mass 
function that can be described with a Schecter mass function 



\ mJ ' 



(37) 



where M c is a cutoff mass around 5 x 10 in starbursts 
ijMeurer et al.l 1 19951 : iMcCradv fc Graham! l2007h . and the 
mass function applies only above a lower mass limit Mi, 
which I take to be 1000 Mq. 
Now suppose that 

• The initial mass function of SSCs has the same form as 
the observed mass function. 

• The ionizing photon luminosity Qi on of a SSC is di- 
rectly proportional to its stellar mass. Furthermore, Qi on is 
assumed to be a reverse step function: a constant for SSCs 
of an age less than t- lon = 10 Myr, and then zero afterwards. 

• The starburst has been continuously forming stars (and 
SSCs) at a constant rate for a time fburst > feon- 

Then I can construct a mass function of SSCs with ion- 
izing stars: 

riiVio, 
dA'h 



= C ion AC 2 exp(-M,/M c ), 



(38) 



which is normalized so that Mi° n — SFR x ti 0n = 
CionM- 1 exp(-M*/M c )dM*. For Mi = 1000 M and 
M c = 5 x 10 6 M , 



Go 



0.126 SFR h 



(39) 



Under my assumptions, these ionizing photons all come 
from stars with ages less than ti on . We can therefore convert 
the star-formation into the mass of the SSCs containing ion- 
izing stars as SFR = Mi° n /t ion : 



Qio 



2.18 x 10 46 ph s" 1 



A I 1 : 



Mr. 



© 



(40) 



Plugging in typical values for an SSC in a starburst, I find 



Rs = 8.8 pc 



AL 



10 6 Mq 



1/3 



II II 



1000 cm- 



-2/3 



(41) 



The absorption coefficient at low v from H II regions is, 
after integrating the SSC mass function from Mi to infinity, 



3Cion7t^0 



VM, 



2/3 



exp(-Af,/M c ) r(2/3, Mi/M c ) 



M, 



1/3 



Al 



1/3 



(42) 



where I take R s = R (M/M ) 1/3 . For the case when Mi 
1000 Mq and M c = 5 x 10 6 Mq, I find 



Qcff 



= (4600 pc)" 



xf UH -) 
V1000 cm" 3 / 



SFR 

Mq yr- 

■4/3 / 



/ Rsb " 
250 pc 



^ion 

10 Myr / 

2 / 



'SB 



100 pc 



|(43) 



The absorption coefficient is smaller than in the case of each 
O star having its own Stromgren sphere. This can be seen if 
we assume each cluster has JV* O stars, each with the same 
ionizing luminosity. Since the Stromgren radius goes only as 



1/3 

2/3' 



the cross section of each H II region increases only as 
The number of Stromgren spheres instead decreases 



as N+ , meaning the effective absorption coefficient is pro- 
r /3 

portional to N+ . clustered stars are not as effective at 
obscuration. Essentially, clustering preserves the filling fac- 
tor that is ionized, but since the regions that are ionized 
are spatially correlated - the ionized regions at the front 
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0.001 



100 




100 



Figure 3. Effective absorption coefficient for H II regions of 
various densities, as normalized to a volumetric star-formation 
rate of 1 Mq yr — 1 kpc — 3 . This plot shows the case when O 
stars are all within 'simple' SSCs surrounded by H II regions, 
in which the ionizing photon luminosity of a SSC remains con- 



stant for t-,. 



10 Myr and then shuts off. The low mass cut- 



off is Mi = 1000 Mq, and the characteristic highest mass is 
M c = 5 x 10 6 M . 

of a large H II region already obscures the ionized volume 
behind it - the covering factor is decreased. 

I plot in Figure [3] the effective absorption coefficients 
for Stromgren spheres around SSCs with masses above 
1000 M and a cutoff of 5 x 10 6 Mq. The ot eS are indeed 
lower than in the case of isolated O stars, but the functional 
form is basically the same, with a plateau at low frequen- 
cies that depends on hydrogen density and an asymptotic 
form at high frequencies that is independent of hydrogen 
frequency. 



3.3 H II regions around Super Star Clusters: A 
model that includes aging 

My assumption in the previous subsection - that SSCs be- 
have like light bulbs, emitting ionizing photons for some time 
before shutting off - is simplistic. In fact, with models such 
as Starburst99, it is possible to predict how the ionizing pho- 
ton g eneration rate evolves for stellar populations of various 
ages l|Leitherer et al.ll 19991 ). In principle, one can then use a 
known star-formation history to accurately predict the clus- 
ter mass and age distribution function, and then integrate 
the cross sections to get an effective absorption coefficient. 

As an example, I will consider the case when the star- 
formation rate has been constant for a duration tburst, before 
which it was zero. The SSC distribution function is then 
assumed to have the form 



dN 



dM+dt i bu 



-M~ 2 exp 



V MJ 



Figure 4. Effective absorption coefficient for H II regions around 
aging SSCs of various densities, as normalized to a volumetric 
star-formation rate of 1 Mq yr — 1 kpc - 3 . The starburst age is 
assumed to be 10 Myr. The low mass cutoff is iUj = 1000 Mq, 
and the characteristic highest mass is M c = 5 X 10 6 Mq. 



where M* is the mass of stars initially formed in the 
cluster, and t is the age of the cluster. The SSC ini- 
tial mass function then has the same form as before, in 
equation 1371 The normalization Caging is again set by 
integrating over masses to get the star-formation rate, 
SFR = Mi<dN/(dM*dt)dM*. For my standard values 

of Mi = 1000 M and M c = 5 x 10 6 M , I find C ag ing = 
0.126 SFRt b 

urst • 

I find the effective absorption coefficient by integrating 
the H II region absorption cross section over different cluster 
masses and ages: 



CKefi = 



dN 1 
dM+dt V 



a{Q ioa {M i ,,t) : u)dtdM t ,. (45) 



The starburst volume is here denoted as V . The ionizing 
photon rate Qi on for a stellar population of age t is given in 
iLeitherer et~al] (|l999l ). 

I show the resulting a c a for a 10 Myr old continuously 
forming starburst in Figure [4] The results are very similar 
to the simply-modelled SSCs (as plotted in Figure [3}. 



4 FITS TO M82'S RADIO SPECTRUM 

The brightest starburst in the radio sky is M82. Lo- 
cated ~ 3.6 Mpc aw ay (as a dopted in this section from 
iFreedman et al.ll 1994 although ISakai fc Madore3ll999l mea- 
sure a distance of 3.9 Mpc), it has a total infrared lumi- 



(44) 



nosity of 5.9 x 10 10 Lq jsanders et al.ll2003T ). correspond- 
ing to a Salpeter IMF star-formation rate of 10 Mq yr _1 
|Kennicuttlll998l ). Most of the radio and infrared emission 
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comes from a region of radius ~ 250 pc (|Goetz et al.lll990l ; 
I Williams fc Bowerll2010h . The starburst is viewed essentially 
edg e-on from Earth. 

IWilliams fe Bowerl (|2010h obtained high quality ra- 
dio observations in the GHz range using the Allen Tele- 
scope Array (ATA). They also compile interferometric and 
single-dish observations in the frequency range of 20 - 100 
GHz, which are useful in constraining free-free emission. 
Noting a systemat i c offs et in the single-dish observations, 
IWilliams fc Bower! (|2010h do not include single-dish obser- 
vations in their modelling, a practice I follow for observa- 
tions above 1GHz. Unlike the other bright starburst, NGC 
253 l|Carillilll996l l. most of the radio emission comes from 
the s tarburst itself rather than the host galaxy (|Basu et alj 
|2012| ). This is especially important at low frequencies, where 
galaxies are frequently unresolved: the starburst may be ob- 
scured by its own H II regions, but the host galaxy is likely 
to be unobscured. 

Besides these observations, there are quite a few ob- 
servations below 1 GHz, although the errors are naturally 
larger (see Table QJ. In theory, the wide frequency cover- 
age of the data, spanning from 22.5 MHz to 92 GHz, makes 
M82 a good choice for spectral modelling. However, I note 
that many different instruments were used in collecting this 
data, with widely disparate beam sizes. For most of the mea- 
surements below 1 GHz, not even the host galaxy (diameter 
~ 10!) is resolved. In contrast, the ATA has a synthesized 
beam diameter of 4.' 2 at 1 G Hz and 35" at 7 GHz, suf ficient 
to resolve the host galaxy ( Williams fc Bowerl [2010:). The 
VLA and GMRT low frequency observations also had beam 
sizes small enough to resolve the host galaxy (and the star- 
burst itself for the GMRT; ICohen et all 120071 : iBasu et all 
l2012f ). The radio flux at 333 MHz from the GMRT is be- 
tween those of the unresolved 178 and 750 MHz measure- 
ments, suggesting that little flux is missing. The radio flux 
from the 74 MHz VLA sky survey is only ~ 2/3 of that 
from the unresolved 57.5 and 86 MHz measurements, which 
may mean that some flux is missing (perhaps from the host 
galaxy). However, the error bars are very large, so it is un- 
clear this is the case; furthermore, the 74 MHz VLA sky 
survey did report i ntegrated fluxes even for resolved sources 
l|Cohen et al.ll2007n . In any case, the 74 MHz VLA sky sur- 
vey reports a major axis size for M82 of 66'.'9± 2'.'8 (smaller 
than the beam size), indicating that the majority of M82's 
74 MHz radio emission comes from within 600 pc of its cen- 
tre (|Cohen et al.ll2007h . Thus, it appears the starburst itself 
is emitting at these frequencies, not just the host galaxy. 

For this paper, I largely ignore the different beam sizes 
and assume all of the radio data points accurately measure 
the radio flux from the inner starbursting region of M82. 
The results should therefore be treated cautiously, until this 
assumption can be verified by high resolution studies at low 
frequency, such as with the GMRT or the Jansky VLA. My 
primary purpose here is to demonstrate how low frequency 
radio spectra can be fit with the new free-free absorption 
models, and to show what kinds of information can be ex- 
tracted. 

4.1 Procedure 

I fit the properties of the H II regions. The star- 
formation rate sets the number of H II regions and is 



allowed to be SFR = 0.625,1.25,2.5,5,10,20 M yr. 
The electron density is allowed to be n e = 
100, 200, 300, 600, 1200, 1800 cm" 3 , and the electron temper- 
ature can be T e = 5000, 7500, 10000, 12500, 15000, 20000 K; 
these set the Stromgren radius and turnover frequency for 
each H II region. Finally, I set the scale height of the star- 
burst hsB, which affects the covering fraction, to be 100 pc. 
I consider both a scenario where the H II regions surround 
individual O stars and where the H II regions surround 
SSCs with a low mass cutoff at 1000 Mq and a high mass 
cutoff at 5 x 10 6 M ; the SSCs are modeled either as simple 
or aging continuously. In the agin g SSC model, I assume 
M82' s starburst is 15 Myr old IjForster Schreiber et al.l 
l2003h . 

I then suppose the radio emission to be a combination 
of synchrotron and free-free emission. Rather than running 
models of cosmic ray e populations, a process which would 
introduce many fre e parameters, I use the p henomenologi- 
cal function form of Will iams fc Bowerl (|2010h for the unab- 
sorbed synchrotron spectrum: 



logi 



cm 



Jy 



= A+B log 



10 (ghz 



1<H2 



GHz 



(4 



In this formula, — B corresponds to the spectral index, 
thought to be near 0.6 — 0.8, and C is the spectral cur- 
vature, thought to be < due to increasing synchrotron 
and Inverse Compton l osses for electrons and high energies 
^Thompson et al.ll2006l ). 

Any H II regions that contribute to free-free absorption 
must also emit free- free emission. The ionizing photon lumi- 
nosity that supports these H II regions, as parametrized by 
the star-formation rate, is directly related to their thermal 



emission as 



niinabs _ 
S mm-ff — 



V n 



I GHz 



= Q22J-,Y SFR V Te )° M ( ° ) 
y \M Q yr- 1 J V 10 4 KJ ^3.6 Mpcy/ 

using the Qi on from equation 1301 [1 Condon! 1 1992h . However, 
it is possible that there is ionized gas which does contribute 
to the free-free emission but not to absorption. This can 
happen if there is a population of small but dense H II re- 
gions with low covering fraction, such as ultracompact H II 
regions. I therefore fit an additional component of free-free 
emission: 

SESfc = ©add (^) Jy. (49) 
The total predicted radio spectrum of M82 is then 



sr d = [5, 



unabs . nunabs 
+ <Jmin-ff 



nunabs 
' "5add-ffJ 



F 



(50) 



Minabs 

I use the edge-on disc value of F/F una ,bs as given in equa- 
tion [23] 

For each combination of parameters describing radio 
absorption (SFR, n e , T e , Hsb), I use \ 2 fitting to find the 
values of A, B, C, and 2\dd that best fit the radio spectrum. 
The values of A are allowed to range from to 2, with a 
spacing of 0.01; I try values of B from -0.9 to -0.3, with a 
spacing of 0.02; C ranges from -0.2 to 0.1, with a spacing of 
0.02; and ©add ranges from to 1 with a spacing of 0.025. 
Then, I compare the \ 2 values for the fits for each absorption 
parameter set to find the best-fitting model. 
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Table 1. Low Frequency Radio Data for M82 



Reference 



Synthesized beam size (8) a Instrument 



(MHz) 


(Jy) 


22.5 


39 ±5 


38 


23 ±3 


57.5 


29 ±6 


74 


18.36 ± 1.86 


86 


27.8 ±3.8 


151 


16.82 c ±0.62 


178 


15.3 ±0.7 


333 


14 ±1 


750 


10.7 ±0.5 



iRoeer. Costain. fc Steward lll986h 66.' x 102.' 

iKellermann. Paulinv-Toth~ fc Williams 45.' X 45'. secf 



llsrael fc Mahonevl dl990l '> 
I Cohen et all <2007D 



Artvukh et al 



7'.0 X 6.'5 



80'.' 



<196Sh : ~ 12.' X 12 .' 



Laing fc Peacock! jl980h 
Baldwin etHTl dl985l) ; lHales et all 4^2 x 4.'2 cosec <5 
jl99ll) 



IKellermann et al.l dl969h 

i Basu et all J2012h 
IKellermann et al.l (|l969f) 



23.' x 18.' sec f 

22'.' x 15'.' 
18.'5 x 18.'5 



Dipole array, Dominion Radio Astro- 
physical Observatory 
Aperture-synthesis system, Mullard 
Radio Astronomy Observatory 
Aperture-synthesis system, Clark Lake 
Radio Observatory 
VLA 

DKR-1000 

6C aperture-synthesis telescope 

Aperture-synthesis system, Mullard 
Radio Astronomy Observatory 
GMRT 

Green Bank Telescope 



a : The zenith angle i s ( and the declina tion is S. 
b : As recalibrated bv lKlein etafl (Il988h . 
c : Average of two measurements. 



As a counterpoint to these models, I also considered 
models where there was no free-free absorption, and uniform 
slab models. In each case, A, B, and C were free parameters. 
For the no-absorption models, I simply added a free-free 
emission component parametrized by T> a dd- In the uniform 
slab models, the electron temperature T e was a free parame- 
ter, with the same values as in the discrete H II region mod- 
els. I then chose T> to fit the free-free emission and absorp- 
tion, where the absorption coefficient is directly related to 
the emission as a v = j v /B v (T e ) = S V D 2 /(27rB„(T e )fl§ B fc), 
using the blackbody source function B v (T e ). In each case, 
I used x 2 fitting to the radio data to select the best-fitting 
parameters. While I fit to all of the radio data for the no- 
absorption models, for the uniform slab case I only used data 
points above 1 GHz. When I include the low frequency data, 
the best-fitting models invariably have no free-free absorp- 
tion or emission. 

I also redid the fits using only the radio data points 
where the beam size was less than 10.', to see how measure- 
ments with large beam sizes were affecting my results. 



4.2 Results 

The best-fitting results with discrete H II regions, whether 
surrounding individual O stars or SSCs, are better at re- 
producing the low frequency radio spectrum of M82 than 
either a fit without absorption or a uniform slab model (Ta- 
ble [2J. The best-fitting individual O star model is shown 
in Figure [S] the best-fitting simple SSC model is shown in 
Figure [6] and the best-fitting aging SSC model is shown in 
Figure [7] While the simple and aging SSC models both have 
low temperatures and densities, the best-fitting individual 
O star model has higher density and temperature. However, 
there are several conclusions that are similar in each of the 
three cases: 

• The unabsorbed synchrotron spectrum is fairly flat, 
with an intrinsic nonthermal spectral index of ~ 0.6, which 



is lower than the t ypical values of ~ 0.8 — 0.9 for normal 
spiral galaxies (e.g.. iNiklas et aj]|l997t ). Previously, this has 
been attributed to free-free absorption, but my models show 
this is not necessarily the case: the spectrum really could be 
intrin sically flat. My values are basically consistent with the 
fit of I Williams fc Bower! (|2010l ), B = -0.56 ± 0.02. 

• The unabsorbed synchrotron spectrum is curved, be- 
coming steeper at higher frequency. The derived intrinsic 
nonthermal spectral curvature in the best-fitting models is 
consistently C ~ —0.08. Curvature is expected if ionization 
and brems strahlung losses are responsible for fla ttening the 
spectrum l|Hummellll99ll : [Thompson et al.ll2006l). This cur - 
vature is slightly less than that of Williams fc Bower! l|2010l ). 
C = -0.12 ±0.03. 

• The covering fraction of the H II regions is of order 
~ 1/2, and the filling fraction ranges from 0.2 to 3%. As 
a result, most of the radio flux is transmitted even at low 
frequencies, with transmission fractions of ~ 60 — 75%. At 1 
GHz, free-free absorption reduces the observed flux by only 
a few percent. 

• These models require low ionizing photon luminosi- 
ties relative to M82's star-formation rate of 10 Mq yr -1 . 
Partly, this is because the low frequency radio spectrum 
is consistent with only small amounts of free-free absorp- 
tion. The ionizing photon luminosity powering the absorb- 
ing H II regions is equivalent to a star-formation rate of 
only < 1 Mq yr _1 . Even when I include an additional non- 
absorbing free-free emission component, the equivalent star- 
formation rate is only ~ 2 Mq yr _1 . The thermal fraction 
at 1 GHz in these models is only ~ 5%. 

In contrast to the discrete H II region models, the best- 
fitting model with no free-free absorption (solid grey line 
in Figures [5] - [7]) requires strong intrinsic spectral curvature 
(C = —0.18) and a flatter synchrotron spectrum (B = —0.52) 
in order to not overproduce the observed low frequency emis- 
sion. However, at the lowest frequencies, the spectral curva- 
ture then suppresses the radio emission too much, leading to 
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Table 2. M82 Radio Spectrum Fits 



Quantity 


No Absorption 


Uniform Slab 


Individual O Stars 


Simple 


SSCs 


Aging 


SSCs 


Data 


A1P 


> GHz 


All a 


All 


e < 10'. b 


All 


< 10.' b 


SFR abs (M Q yr~ 1 ) c 
n e (cm -3 ) 






0.625 


0.625 


0.625 


1.25 


1.25 




35 


600 


100 


300 


100 


300 


Te (K) 




20000 


20000 


7500 


20000 


7500 


20000 


ft (pc) 


100 


100 


100 


100 


100 


100 


100 


.4 


0.91 


0.95 


0.93 


0.94 


0.92 


0.94 


0.92 


8 


-0.52 


-0.60 


-0.58 


-0.60 


-0.56 


-0.58 


-0.56 


C 


-0.18 


-0.06 


-0.10 


-0.08 


-0.14 


-0.08 


-0.14 








0.19 


0.12 


0.19 


0.23 


0.36 


23 add 6 


0.65 


0.30 


0.33 


0.35 


0.45 


0.13 


0.28 


x 2 


96.0 (75.5°) 


69.7 (657^) 


83.0 (75. 8 a ) 


81.8 


74.6 


81.6 


74.1 


/cover 




100% 


43% 


50% 


25% 


57% 


29% 


/fill 




100% 


0.23% 


2.3% 


0.57% 


2.8% 


0.70% 


/ thcrm (1 GHz)/ 


7.4% 


3.4% 


6.1% 


5.4% 


7.2% 


4.1% 


7.1% 


^/-FunaUf -+ 0) 9 




0% 


74.6% 


61.9% 


81.7% 


64.8% 


83.4% 


^/FunabsCl GHz) h 




96.7% 


97.8% 


96.5% 


97.9% 


97.3% 


98.4% 


Pan/k (K cm" 3 ) 1 




6.9 x 10 5 


3.6 x 10 7 


2.3 x 10 6 


1.8 x 10 7 


2.3 x 10 6 


1.8 x 10 7 


SFRem (M yr" 1 ) 


2.9 


0.99 


1.7 


2.4 


2.1 


1.8 


2.1 



a : The best-fitting parameters for these models are the same whether all data are included or just those with beam sizes less than 
10 arcminutes. The x 2 value when only data with beam sizes less than 10 arcminutes is given in parentheses. 
b : Fluxes for observations with beam sizes less than 10 arcminutes (see Table [TJ. 

c : Effective star-formation rate, which sets the ionizing photon luminosity of the starburst. This sets both the number of H II 
regions and the amount of free-free emission. 

d : Thermal free-free emission at 1 GHz from the H II regions responsible for free-free absorption. This sets a floor on the amount 
of free-free emission. 

e : Additional free-free emission from ionized gas that does not contribute to free-free absorption in my model (e.g., from more 

compact, denser H II regions). This is a free parameter. 

': Fraction of the flux at 1 GHz which is thermal free-free emission. 

9 : Fraction of the radio flux transmitted by the free-free absorbing H II regions at extremely low frequencies. The effects of free-free 

absorption in other phases or synchrotron self-absorption are not included. 

h : Fraction of the radio flux transmitted by the free-free absorbing H II regions at 1 GHz. 

': Derived thermal pressure in H II regions, 3n e T e . 

i : In selecting the uniform slab model, x 2 was minimized for data above 1 GHz; this value is 69.7. However, if all radio data is 
included, the total \ 2 increases to 657. 



a relatively poor fit. The radio flux in the uniform slab model 
(long-dashed line) plummets rapidly and cannot explain the 
data points below ~ 300 MHz at all. 

The effects of free-free absorption on the spectrum 
shape are concentrated within a finite frequency band. At 
high frequencies, the H II regions are transparent. At low fre- 
quencies, the H II regions are completely opaque, and a e g 
reaches some constant value; this changes the normalization 
of the spectrum, but not its intrinsic shape. These effects 
can be described more quantitatively in terms of the total 
spectral index B and spectral curvature C: 



B 



C 



dio gl0 sr d 

d log 10 v 

d 2 i ogl0 sr d 



(51) 



(52) 



d ( lo Sio^) 2 

these quantities differ from B and C in that they include 
free- free emission and absorption. I show how B and C vary 
with frequency in Figure [H] for the best-fitting modelfl Free- 
free absorption (the difference between black and grey lines) 
flattens the spectrum (higher B) between ~ 0.1 — 1 GHz: this 



5 Note the use of base 10 logarithms. While B = d In Sj rcd 
C = In 10 X d 2 In Sg rcd /d(ln u) 2 



/dlnu, 



is the regime where H II regions become optically thick. The 
resulting pulse in B is narrower when H II regions surround 
individual O stars (dotted lines) instead of SSCs: in that 
model, the H II regions all have the same size, with the 
same turnover frequency, instead of the range of turnover 
frequencies for SSCs. A pulse also appears in the spectral 
curvature in this frequency range. One potential issue for 
studies looking for intrinsic nonthermal curvature in the ra- 
dio spectrum is that either free-free absorption or emission 
affect most observable radio frequencies. Only at frequencies 
of a few GHz or, paradoxically, below ~ 30 MHz where cv e ff 
is constant, does the intrinsic curvature dominate over that 
introduced by free-free absorption or emission. 

One motivation for studying free-free absorption is to 
calculate the thermal press ure in H II re gions from the den- 
sity and temperature (c.f.. ICarillilll99f3 ). Even in my best- 
fitting models, the pressure varies by an order of magni- 
tude. It is extremely high in the individual O star model, 
with P/k = 4 x 10 7 K cm . In my SSC models, it is only 
P/k = 2 x 10 6 K cm . Numerous studies have inferred the 
pressure of M82's H II re g ions, using infrared and optical 
spectroscopy. ISmith et al.l (|2006l ) found a fairly high pres- 
sure of P/k = (l-2)xl0 7 K cm for the H II region around 
the M82 A-l super star cluster, still less than the individ- 
ual star model, but much greater than in the SSC models. 
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Figure 5. The effects of free-free absorption on M82's radio spec- 
trum, for H II regions surrounding individual O stars. The solid 
black line is the best-fitting to the free-free absorption, and the 
dotted line is what that fit spectrum would be without absorp- 
tion. The grey solid line is the best possible fit with no absorption, 
though it does include free-free emission. I show a uniform slab fit 
to the > f GHz radio spectrum as the dashed line. The residuals 
are plotted on the bottom. 
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Figure 6. A model of the free-free absorption on M82's radio 
spectrum, using the simple super star cluster assumptions de- 
scribed in section 13.21 I plot the best-fitting model to the radio 
data with beam sizes < 10' as the dash-dotted line. The other 
line styles are the same as in Figure [5] 
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Figure 7. A model of the free-free absorption on M82's radio 
spectrum, where the model now accounts for how the ionizing 
photon luminosity of SSCs evolve as they age. I plot the best- 
fitting model to the radio data with beam sizes < 10' as the 
dash-dotted line. The other line styles are the same as in Figure[5] 



Note that M82 A-l is relatively large (~ 10 6 Mq) compared 
to the small SSCs t hat would be expected to d ominate the 
free- free absorption. IWestmoquette et al.l l|2007l ) found lower 
pressures of P/k = (5 - 10) x 10 6 K cm" 3 in most M82 H II 
regions, though these pressures are st ill several times la rger 
than what I find in my SSC models. lLord et"all (| 1996T ) in- 
ferred H II region pressures of P/k = 3 x 10 6 K cm -3 , based 
on infrared diagnostics of the surrounding denser and colder 
photodissociation regions. In any case, though, because of 
the range of densities I allowed for discrete H II regions, the 
pressures are necessarily higher than what I would find with 
a uniform slab model (~ 7 x 1 5 K cm" 3 ; c ompare with the 
similar results for NGC 253 in lCarilli|[l99rj) . 

The low ionizing photon luminosities I derive are wor- 
rying, especially since I allow for 'hidden' free-free emission 
from compact H II regions that do not contribute to ab- 
sorption. The low luminosities may be because of the large 
errors and sparseness of the high frequency data. On the 
other hand, thermal dust emission probably contributes to 
the highest frequency data points at some level, and there 
may be spinning dust emission as well, but including any 
dust emission would just tighten the constraints on the free- 
free emission. It is possible, though, that the amount of free- 
free emission really is low in starburst galaxies - perhaps 
not a surpr ising hypothesis given that st arbursts are dusty 
places (c.f.. |Petrosian. Silk, fc Fieldlll972T ). Or perhaps some 
ultraviolet photons escape into the hot superwind phase that 
likely fills much of the starburst volume, instead of ionizing 
the molecular gas. If so, then free-free emission is not neces- 
sarily a simple star-formation rate indicator in starbursts: its 
strength depends somehow on the radiative transfer in dusty 
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its high sensitivity, the shape of the spectrum can be more 
meaningfully constrained. Moreover, as a single instrument, 
there would be fewer worries about systematics when com- 
paring between data at different frequencies. Furthermore, 
it has high spatial resolution, so it can actually measure the 
characteristic sizes of H II regions, in particular searching 
for small, high density H II regions. Another possibility is 
to use radio recombination lines , also observable with the 
Jansky VLA l|Keplev et al.|[201ll) . which constrain not only 
the ionizing photon luminosity of a starburst, but which 
density and temperatures t he ionized gas has. I note that 
iRodriguez-Rico et all ([2004) found using H53a and H92a 
radio recombination lines an ionizing photon luminosity of 
1.6 x 10 53 s _1 for M82, equivalent to a star-formation rate 
of 0.7 Mq yr _1 , compatible w ith the low ionizing photons I 
find here. On the other hand, iPuxlev et al.l (|l989T ) deduced 
an ionizing luminosity of 1.1 x 10 54 s _1 also using the H53a 
line, which would correspond to a star-formation rate of 
5 Mq yr _1 . This is significantly greater than what I find 
with my fits, but still a factor of ~ 2 lower than the infrared- 
derived star-formation rate. Future radio recombination line 
studies may clarify the situation. 



Figure 8. How free-free absorption affects the total spectral in- 
dex B (top) and total spectral curvature C (bottom). I show these 
values with (black) and without (grey) free-free absorption. The 
line styles are for different M82 models; dotted is individual O 
stars, dashed is simple SSCs, and solid is aging SSCs. Free-free 
absorption introduces a pulse in the spectral index and curvature, 
as H II regions transition from being transparent to opaque with 
decreasing frequency. 



environments, and it is buried even more than expected by 
synchrotron emission. 

As I noted before, though, these results must be con- 
sidered preliminary, since the low frequency data so vital 
in constraining free-free absorption comes from many in- 
struments with different beam sizes. A future low frequency 
survey of M82 that resolves the starburst would ensure that 
flux from the surrounding galaxy is not making M82's star- 
burst appear brighter than it really is; such a survey can be 
done with the GMRT, VLA, or LOFAR. 

When I fit only the data with beam sizes smaller than 
10.', the basic conclusion that most of the flux is transmit- 
ted stands. The best-fit individual O stars model is actu- 
ally unchanged with these data. However, the models of 
Stromgren spheres around SSCs now fit better for somewhat 
higher densities (300 cm -3 ) and much higher temperatures 
(20000 K). The best-fitting spectral index is now 0.56, with a 
curvature C = —0.14, higher than before. These models have 
even lower H II region covering fractions (25 — 30%), but 
comparable thermal fractions (5 — 7%). Therefore, the en- 
vironments of the H II regions are not too well constrained, 
and better radio data is needed, but absorption from dis- 
crete H II regions are adequate so far to explain the radio 
spectrum. 

Besides better low frequency data, it would be helpful 
to have high quality observations at 10 - 100 GHz. Then it 
may be possible to meaningfully detect the flattening from 
free- free emission, if it is pres ent. The Jansky V LA is able 
to observe at up to 50 GHz (|Perlev et al.ll2011u , and with 



5 WHERE DO THE RADIO SPECTRA OF 
STARBURSTS REALLY END? 

The premise of this work is that free-free absorption from H 
II regions only partially covers starburst galaxies, and star- 
bursts still emit quite a lot of radio emission at low frequen- 
cies. However, the starburst radio spectrum cannot continue 
indefinitely to zero frequency. The presence of more diffuse 
volume-filling phases of ionized gas within starbursts will 
provide free-free absorption. In addition, the purely nonther- 
mal process of synchrotron self-absorption eventually must 
cut off the radio spectrum of starburst galaxies even if there 
is no diffuse ionized gas. In this section, I calculate the ef- 
fects of these additional processes to find when starbursts 
should have true cutoffs in their radio spectra. 

5.1 Free- free absorption in starburst winds 

The high rate of supernovae in starburst galaxies occur- 
ring within a relatively small s pace is expected to ex- 
cavate a hot phase of the ISM (|McKee fc Ostrikerl Il977l ; 
iHeckman. Armus. fc Milevl [l990h . The hot ISM erupts as 
a starburst-wide superwind, and is low density, but high 
in pressu re, temperature, and, it is thought, filling fac- 
tor (e.g., IChevalier fc Cl cgg 1985). Gas emitting in soft 
X-rays is indeed observed in many starbursts, though 
IStrickland fc Stevens! (|2000T ) argues that this emission comes 
from a cooler phase with lower filling factor than the actual 
wind. In addition, Chandra detected diffuse 6.7 keV iron 
line emission in M82 that suppo rts the presence of 10 8 K 
gas (Strickland & Heckman 2007). 

IStrickland fc Heckmanl (]2009) give the central density 
of the superwind as 



^3/2 

Pc = 0.930^ 



M3/2 



therm v ' 



(53) 



and the central temperature for a completely ionized wind 
as 
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T c = 



0.8m H Cthcrm E 



(54) 



k/3M 

In these equations, e t herm ~ 0.5 is the fraction of supernova 
energy that ends up in therm al energy of the plasma, j3 ~ 2 
is the mass- loading fractiorjf] (Strickland & Hcckman 2009). 
They give the mass injection rate from supernovae and stel- 
lar winds as 



M = 0.14 M yr" 1 



SFR 



(55) 



and the energy injection rate from supernovae as 

E = 2.9 x 10 41 erg s" 1 ( ^ SFR \ (56) 

where I have converted between SFR above 1 Mq and the 
SFR for a Salpeter IMF from 0.1 to 100 M using SFR(> 
1 Mq) = 0.45 SFR. Then, in numerical terms, I have 



kinetic energy to radiative losses as they expand in a 
dense molecular medium, so that supernova-heated mate- 
rial fills only small isolated bubbles. Then the molecular 
gas will fill most of the starburst region, since the H II 
regions also have a small filling factor. Molecular gas is 
mostly neutral and dust extinction will rapidly extinguish 
any ultraviolet light. However, cosmic rays should provide 
a relativ ely high level of ionization through fairly large 
colum ns (|Suchkov. Allen, fc Heckmanlll993l ; Papadopoulosl 
2010), though the details of CR diffusion in starbursts and 
their penetration into molecular material is poorly under- 
stood. 

As I argued in lLackil (|2012h . the cosmic ray ionization 
rate in starburst galaxies, where the proton spectrum is rel- 
atively hard, is 



Ccr 



n c = 0.011 cm 
and 



-3 (§_ 

2 



3/2 ft^ \ - 1 ' 2 

c tnerm * 



0.5 



Mq yr -1 kpc 



T c = 8.3 x 10 7 K (^^y 1 ) (|) 



(58) 



I find that the high temperatures and low densities of 
starburst winds makes them extremely poor at free-free ab- 
sorption. The free-free optical depth is tiny even at MHz 
frequencies: 



4.2 x 10" 



-'SFR 



50 M Q yr" 1 kpc 



(l MHz) 



(— 

\ 100 pc 



(59) 



assuming that gg = 20 for MHz frequencies and T m 10 8 K 
plasmas. The turnover frequency is far below observability, 
since we are stuck in the Milky Way: 



6.5 kHz 



50 Mq yr" 1 kpc -2 / V 100 P c 



1/2 



(60) 



The weakest starbursts have Esfr ~ 1 Mq yr 1 kpc 2 
and the star-formation rate surface density in M82 and 
NGC 253 is ~ 50 M yr -1 kpc -2 , while Ssfr reaches 
~ 10 4 Mq yr -1 kpc" 2 in the most extreme starbursts. 
Therefore, the starburst wind introduces free-free absorp- 
tion only at frequencies less than ~ 0.1 — 1000 kHz, and is 
not important at the observable MHz radio frequencies fit 
by models. 



5.2 Free- free absorption from cosmic ray-ionized 
molecular gas 



iThompson. Quataert. fc Murray! l|2005l ) have argued that 
cold molecular gas, instead of rarefied supernova-heated 
material, fills most of the volume of starbursts. This 
could happen if supernova remnants rapidly lose their 

6 The fraction of mass ejected by stellar winds and supernova 
that ends up in the wind; it can be greater than 1 if cold gas is 
swept up by the wind as it leaves the galaxy. 



'/ion 

E ion M H 



(61) 



In this equation, Lcr is the luminosity of injected cosmic 
(57) rays, £"ion ~ 30 eV is the energy lost p er cosmic ray ioniza- 
tion event (jCravens fc Dalgarndll978r i. Mh is the mass of 
gas in the galaxy, and 77i on ~ 0.1 is the fraction of cosmic 
ray power that goes into ionization. The injected cosmic ray 
luminosity is though to be approximately 10 50 erg per su- 
pernova. For a Salpeter IMF extending from 0.1 to 100 Mq, 
the supernova rate is r SN = 0.0064 yr _1 (SFR/M Q yr" 1 ). 
Therefore, the cosmic ray luminosity is approximately 
Lcr w 2.0 x 10 40 erg s -1 (SFR/M yr" 1 ). 

Note that £cr oc SFR/Aftr, the inverse of the gas 
consumption time r gas = Mh /SFR. This time is equal to 
~ 20 Myr for nuclear starbursts, so the ionization rate can 
be conveniently expressed in terms of r gas . The cosmic ray 
ionization rate is: 



Ccr w 1. 



x HT 15 s" 1 



20 Myr 



(62) 



The ionization fraction is set by the ratio of the 
gas density nn and a characteristic density n c h ~ 
1000 (Ccr/10~ 17 s" 1 ): 



x e = lxl0- 7 (^ 



1/2 



V ins J KAuhJ 



(63) 



from iMcKeeJ (|l989l ). In the cosmic ray ionized gas of star- 
bursts, 7i c h 2> tih and x e as 10" 



3.5 x 10" 



20 Myr 



)"( 



7 {n ch /n H )- Thus, 
nn 



500 cm- 



-i 



(64) 



While the ionization fraction is low, free-free absorp- 
tion is enhanced by two factors: molecular material is dense, 
so the density of electrons and ions is relatively high, and 
molecular gas is cold. The electrons and ions reach ther- 
mal equilibrium with the surround ing gas long before they 
recombine (e.g.. iMcCall et al]|2002l ). Assuming typical star- 
burst molecular gas temperatures of ~ 100 K, I find free- free 
optical depths of 



Tff 



17 



20 Myr 



.100 K 



-3/2 



100 pc 



( — 

I MHz 



(65) 



assuming n e = m = x e nn and gg — 10. The frequency of 
the free-free spectrum turnover is 
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vg « 4.1 MHz 



20 Myr 



(100 k) 



-3/4 



100 pc 



1/2 



(66) 



If there really is a molecular medium with a high filling 
factor, and if it is cosmic ray ionized, it can terminate the 
low frequency radio spectrum of starbursts. With turnover 
frequencies of ~ 5 MHz, detecting it would be a challenging 
measurement, but one that is possible in principle, especially 
since the turnover is at higher frequency than that of the 
warm ionized medium in the Milky Way (and presumably 
the starburst's host galaxy), ~ 3 MHz. Thus, the presence or 
absence of a turnover at ~ 5 MHz can test whether cosmic 
ray ionized high filling-factor molecular gas really does exist 
in starbursts. 

A factor that I have neglected is the highly supersonic 
turbulence present in starburst galaxies, which will generate 
large density fluctuations in the molecular gas, even without 
the existence of any other gas phases. However, if the ioniza- 
tion rate £cr across the starburst is the same everywhere, 
the homogeneous approximation will give approximately the 
right results. This is because x e oc and ami oc x^.n 2 H : the 
free-free absorption coefficient does not vary with density. It 
is possible that the cosmic ray ionization rate decreases in 
high densit y regions b ecause the cosmic rays are stopped 
by the gas (|Lackill2012l ). but these regions are those of low 
filling factor anyway. 



5.3 The Razin effect 

Free-free absorption is not the only process that cuts off 
the radio spectrum at low frequency. At low energies, the 
index of refraction of plasma suppresses the beaming of 
synch rotron radiation and caus es it to fall off exponen- 
tially (|Rvbicki fc Lig htman 19791). The frequency where this 
Razin effect becomes important is 



v R = 185 kHz 



(l cm" 3 ) 



( B 



(67) 



^100 /xG / 

from lSchlickeiserl (j2002T l. 

Suppose all the gas in the starburst, both ionized and 
neutral, has the same magnetic field. The Razin cutoff in 
the hot superwind phase is 



51 kHz I ^ SFR , ) ( — *L 
. n.a ) V 50 M yr- 1 J \200 /iG / 

SsFR \ / B ' 



VR = 



340 kHz 



10 4 M yr" 1 J V6 mG 



(68) 



where the upper values are scaled to values expected in M82 
and the lower values are scaled to values expected in Arp 
220's radio nuclei. In fact, B and Esfr are likely to depend 
on each other. The existence of the linear far-infrared radio 
correlation of galaxies constrains magnetic field strengths to 
be 



Bfrc » 67 ^G x /Ssfr/(Mq yr" 1 kpc" 



(69) 



from lLacki et all ( |2010f >. after using the iKennicuttl (119981 ) 
Schmidt law to convert between gas surface density and 
star-formation surface density. Applying this value of B to 
equation [68] the Razin cutoff is 



vW nd = 21 kHz 



50 M yr- 1 



1/2 



(70) 



varying between 3 and 300 kHz for starbursts with Esfr 
between 1 and 10 4 Mq yr -1 kpc -2 . 

If starbursts are instead filled with cosmic ray ionized 
gas, the low electron density (n e = x e nn) of these regions 
imply still lower Razin cutoffs: 



1.6 kHz , 

moi ) V 20 Myr 
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20 Myr 



200 fiG 
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(e mc) 



(71) 



In high density molecular regions of material, the ion- 
ization rate is lower, and Zeeman splitting measure- 
ments indicate the magnetic field may be even higher 
l|Robishaw. Quataert. fc Heilesl 120081 ). so the Razin effect 
could be even less important. 

Where the Razin effect might be important is in H II 
regions, which are fully ionized and high density: 



HII 



46 MHz 



31 MHz 



( — \ (-*—) 

V500 cm- 3 / ^200 pG J 
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(72) 



.10000 cm- 3 , 

Within these regions, we see that free-free absorption - 
which turns H II regions opaque at GHz frequencies - is 
more important than the Razin effect. 

I conclude that the Razin effect is not observable in 
starburst galaxies. 



5.4 Synchrotron self-absorption 

Both free-free absorption and Razin cutoff depend on the 
density distribution of ionized matter, something which is 
not entirely certain in starburst galaxies. However, because 
we observe synchrotron emission from starburst galaxies, 
there must also be synchrotron self-absorption in them as 
well. 

The maximum brightness temperature of a synchrotron 
source, as limited by synchrotron self-absorption, is 



(ghz) 



1/2 



B 



-1/2 



(73) 



100 fiG 

from lBegelman. Blandford. fc Reesl (|l984l ). 

For resolved starbursts with an observed low frequency 
radio spectrum, it is possible to simply fit a model to the 
radio spectrum and see when, if ever, the brightness tem- 
perature T;, is greater than T max . In general, though we ex- 
pect starbursts to lie on the far-infrared radio correlation, 
with vL v (l GHz) = 1Q- 6 L T ir = 5630 WSFR/ M n yr- 1% 



Kennicuttl [1998 : lYun. Reddv. fc Condoni 1200 ll ). Ignoring 



geometrical factors, the starburst can be thought of as 
a sphere with radius J?sb, so that the intensity is I v — 
L u /(4n 2 Rg B ). Likewise, the star-formation surface density is 
^SFR = SFR/(-7riig B ). Finally, since we are in the Rayleigh- 
Jeans limit, Tb = c 2 I^/(2u 2 k). Putting all these formulas 
together, I find that the FIR-radio correlation implies 



T b (l GHz) w 600 K 



JSFR 



Mq yr -1 kpc 



(74) 
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Now we need to extrapolate this down to low radio fre- 
quencies. The simplest assumption is that the radio spec- 
tral index is constant and equal to —0.7. Then by setting 
T b = T b (l GHz)(!//GHz) -2 - 7 , using eon.l69lfor B, and equat- 
ing eqns. [73] and 1741 I find the synchrotron self-absorption 
turnover for M82 is at: 



z, SSA = 3.5 MHz 



50 M Q yr" 1 kpc" 



B 



200 ^G 



(75) 



Plugging in the B field strength from the far-infrared radio 
correlation feqn. Raj?)). I find 



vssa = 4.0 MHz 



50 Mq yr- 1 kpc" 



(76) 



For starbursts with 1 < E S fr/(M yr -1 kpc~ 2 ) < 10 4 , 
i/ssa ~ 0.9 — 30 MHz. This result is conservative, though, 
since the radio spectra display both partially-covered free- 
free absorption and intrinsic spectral curvature. Moreover, 
starburst radio spectra tend to be somewhat flatter than 
a = 0.7. 

I can use my model fit to the spectrum of M82 to 
account for these effects. Taking B = — a = —0.6 and 
C = —0.1, I combine eqns. [74l and 1731 to find 



, / t'SSA \ 



-15.5 + 5\ 5.69 



-0.41og 10 (/S^ F 4 R ), (77) 



where / is the fraction of flux transmitted by free-free ab- 
sorbing H II regions at low frequencies. When / = 1, I find 
synchrotron self-absorption turnovers at 0.3 — 20 MHz for 

1 < Esfr/(Mq yr -1 kpc -2 ) < 10 4 , with a turnover at 

2 MHz for M82; when / = 0.5, the synchrotron turnovers 
for starbursts are nearly the same (0.2 — 20 MHz; 1.5 MHz 
for M82). 

The synchrotron self-absorption cutoffs are higher than 
the free-free absorption cutoffs for the starburst wind. Since 
they are in the MHz range, they may be observable, espe- 
cially for Arp 220 with a frequency cutoff at ~ 5 — 20 MHz. 
A measured synchrotron turnover frequency would tell us 
what the magnetic field strength in starburst galaxies really 
is. 

The Razin effect actually suppresses not only 
synchrotron emission , bu t synchrotron absorption 
|Crusius fc Schlickeiser] [T988:) . Below the Razin cutoff, 
there ca n even be negative absorp tion and stimulated 
emission l|Crusius fc Schlickeiserlll988r i. However, the Razin 
cutoff is at much lower frequencies than the synchrotron 
self-absorption turnover. Furthermore, when the Razin 
effect does operate, there is still a cutoff in the observed 
synchrotron radio emission, so either way the radio 
spectrum ends. 



5.5 Summary 

Of the processes I considered, the most likely to produce 
a true low frequency turnover in starburst radio spectra is 
synchrotron self-absorption, which should be important only 
below ~ 20 MHz even for Arp 220-like compact starbursts. 
If there is a volume-filling molecular medium that is ion- 
ized by cosmic rays, it may also cause a free-free absorp- 
tion turnover, but only below ~ 5 MHz. The volume-filling 



10 K superwind is too hot and rarefied to produce appre- 
ciable free-free absorption at observable frequencies. Nor is 
the Razin effect important. 

It is in principle possible that there are diffuse ionized 
phases that I do not consider, perhaps like the warm ionized 
medium in the Milky Way. Such phases could be constrained 
by the amount of free-free absorption. In fact, the Milky 
Way's own warm ionized medium will prevent observations 
below a few MHz, and it is likely the warm ionized medium 
in starbursts' surrounding host galaxies will do the same, 
especially if they are viewed edge-on. However, given our 
current understanding of the phases likely to fill most of the 
starbursts, it appears that starbursts are largely transparent 
to < 20 MHz. 



6 CONCLUSIONS 

With the renewed interest in low frequency radio astron- 
omy, it is time for a better theoretical understanding of 
the low frequency radio spectra of starburst galaxies. Pre- 
vious models of the emission at these frequencies, if they 
considered free-free absorption at all, predicted that star- 
bursts are opaque below a GHz. This was because they used 
the uniform slab model, essentially assuming that ionized 
gas evenly pervaded the starburst. Most of the ionized gas 
mass is not truly diffuse in starbursts, but instead resides 
in discrete H II regions. In fact discrete H II regions are re- 
sponsible f or much of the free - free absorption i n the Galac- 
tic Centre l|Brogan et al. 2003; lNord et al.l [2006) , a region in 
some ways anal ogous to starbursts but close enough to be 
easily resolved (|Crocker et al.l [201lh . A low frequency map 
of M82 also indicates discrete H II reg ions as the source 
of free- free absorption (| Wills et al.lll997f ). If the H II re- 
gions are uniformly distributed throughout the starburst, 
they contribute to an effective absorption coefficient which 
should be used in the uniform slab formula. This coefficient 
does not approach infinity at low frequencies. H II regions 
only partially cover the starburst, leaving sightlines where 
the emission is unobscured. Furthermore, the H II regions 
are usually some way in from the surface of the starburst, 
so there is unobscured synchrotron emitting material on a 
sightline in front of the nearest H II region. 

I showed how to calculate this effective absorption coef- 
ficient in the limit that the H II regions do not fill most of the 
starburst, which is usually true. The calculation ultimately 
reduces down to the density of the H II regions multiplied by 
the cross section of each to absorb background radio waves. 
I applied the calculations to the radio spectrum of M82. I 
find that models with discrete H II regions, around either 
individual O stars or Super Star Clusters, are able to repro- 
duce the radio spectrum reasonably well, even at frequencies 
down to 23 MHz (see Figures and [7$. 

The models I presented in the paper were relatively sim- 
ple, assuming that H II regions are characterized by a single 
temperature and density. However, if the temperature and 
density distribution is known through other means, such as 
radio recombination lines, it is straightforward to integrate 
up the cross sections and calculate the effective absorption 
coefficient. A more important flaw in the approach described 
here is that I assume that H II regions are uniformly dense 
Stromgren spheres. However, dust absorption affects both 
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the amo unt of ionizing photo ns available to produce H II 
region s (IPetrosian et alj|l972h . and the structure of H II re- 
gions (|Drainei2011bl ). Since my models suggest a low amount 
of free-free emission, which could mean that dust is absorb- 
ing ionizing photons, it is important that these effects be 
studied. 

This work implies that starbursts are in fact fairly 
bright at low frequencies. I considered whether any other 
absorption process could cause a turnover in the radio spec- 
trum. I found that free-free absorption from the volume- 
filling superwind phase of starbursts, and the Razin effect, 
are negligible down to kHz frequencies in starbursts. Syn- 
chrotron self-absorption should become important at fre- 
quencies of a few MHz. If there is a volume-filling phase 
of cold molecular gas, ionized by cosmic rays, it can also 
cause free-free absorption at a few MHz. 
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